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We study properties of the proton component of neutron star mat- 
ter for a number of realistic nuclear models. Protons which form a few 
percent admixture tend to be localized in potential wells corresponding 
to neutron matter inhomogeneities created by the protons in the neutron 
medium. We calculate the energy of the Wigner-Seitz cell enclosing a single 
localized proton. The neutron background is treated in the Thomas-Fermi 
approximation and the localized proton is described by the Gaussian wave 
function. The neutron density profile is obtained by solving the appropri- 
ate variational equation. This approach gives lower energies of localized 
protons than obtained previously with less sophisticated methods. 

PACS numbers: 21.65.+f,97.60.Jd 

1. Introduction 

Physical properties of the neutron star interior relevant to macroscopic 
observables are rather sensitive to the the microscopic structure of dense 
nuclear matter in neutron stars. For example, the transport and magnetic 
properties of neutron stars depend strongly on the structure of the so called 
liquid core. Particularly important is the structure of the proton compo- 
nent. A uniform proton distribution and a periodic (crystalline) proton 
arrangement result in very different properties [1]. The latter possibility 
was discussed in Refs. [2, 3] for strongly asymmetric nuclear matter which 
was shown to be unstable with respect to proton localization. The local- 
ization effect is a result of the interaction of protons with small density 
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oscillations of the neutron background [4]. The protons behave as localized 
polarons which form a periodic lattice at high densities [5] . 

The presence of the localized protons inside neutron star cores would 
have profound astrophysical consequences. The transport properties of such 
a phase are rather different from that of a uniform nuclear matter [6]. In 
particular, the cooling proceeds in a quite different way. Recent analysis [6] 
shows that the presence of such localized proton phase results in more sat- 
isfactory fits of temperatures of observed neutron stars. Also, spin ordering 
of localized protons could strongly affect magnetic properties of the system 
[2, 7]. The spin ordered phase can contribute significantly to the observed 
magnetic moments of neutron stars [8, 9]. 

The aim of this paper is to study the proton localization for a number 
of realistic nuclear models with improved variational method. In original 
calculations [3, 4] both the proton wave function and the neutron density 
distribution were assumed to be trial functions which included variational 
parameters. In this paper we find better estimates of energies of localized 
protons by solving the appropriate variational equation for the neutron den- 
sity profile that gives the minimum energy for a fixed wave function of the 
localized proton. 

The paper is organized as follows: In the next section we describe the 
model of proton impurities in the neutron star matter. In Sect. 3 simple es- 
timates of the proton localization based on trial functions are given. Sect. 4 
contains the formulation and the solution of the variational problem. Re- 
sults are discussed in Sect. 5. 

2. Model of proton impurities in the neutron star matter 

The amount of protons present in the neutron star matter, which is 
charge neutral and /3-stable, is crucial for the cooling rate of neutron stars 
and also plays an important role for magnetic and transport properties of 
neutron star matter. Nuclear models do not uniquely predict the proton 
fraction of the neutron star matter at high densities. This controversy is 
discussed in details in Ref. [1, 10] where the discrepancy of the proton 
fraction in various models is shown to reflect the uncertainty of the nu- 
clear symmetry energy at high densities. In this paper we consider a class 
of nuclear interaction models for which the proton fraction is of the or- 
der of a few percent and decreases at high densities - as shown in Fig.l. 
For the calculations we have chosen six realistic nuclear interaction mod- 
els. These are interactions derived by Myers and Swiatecki [11] (MS), the 
Skyrme potential with parameters from Ref. [12, 13] (Sk), the Friedman and 
Pandharipande interactions [14] (as parametrized by Ravenhall in Ref. [15]) 
(FPR) and three models, UV14+TNI, AV14+UVII and UV14+UVII, from 



szm printed on February 1, 2008 



3 



Ref.[16] by Wiringa et al. 

Let us consider a neutron star matter containing a small proton fraction 
x. To compare the energy of a normal phase of uniform density and a phase 
with localized protons we apply the Wigner-Seitz approximation and divide 
the system into cells, each of them enclosing a single proton [2, 4]. For 
simplicity, the cells are assumed to be spherical. The volume of the cell is 
V = l/np. The normal phase is of uniform density tin and the neutron 
chemical potential is fijy- In the uniform density phase protons are not 
localized and their wave functions are plane waves. 

The energy of the cell, which is a sum of proton and neutron energies, 
reads 



where e(njv,np) is the energy density of the uniform phase. For small 
proton density, i.e. for low x, we can expand the energy density 



In the following we shall adopt abbreviations e (tin) = e (njv,0) for the 
energy density of pure neutron matter and fip (tin) = fJ-p(riN,0) for the 
proton chemical potential in pure neutron matter. The energy of the cell is 
thus approximately 



Our aim is to compare the energy of the normal phase, where protons 
are nonlocalized, with the energy of a phase where the protons are trapped 
into potential wells, corresponding to the nonuniform neutron density distri- 
bution, which most likely form a regular arrangement. We treat this proton 
"crystal" in the Wigner-Seitz approximation. 

Let us consider a Wigner-Seitz cell with nonuniform neutron matter 
distribution n (r) surrounding the proton whose wave function is \t p. In the 
local density approximation one can identify the proton effective potential 
with the local proton chemical potential np{n) [2]. The proton's effective 
potential varies locally with neutron matter density n (r). This results in a 
potential well fxp (n (r)) which affects the single proton wave function. The 
energy of the Wigner-Seitz cell, El, is: 



E = Ve (n N ,n P ) , 



(1) 



e (njv, n P ) e (n N , 0) + fi P (n N , 0) n P . 
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E = iip (n N ) + Ve (n N ) . 
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The first term is the energy of the proton in the effective potential 
v e ff(r) = fip(n(r)). It is by construction the attractive potential well. 
At high densities the derivative of the proton chemical potential is positive, 
> 0, for all interactions we use. This can be seen in Fig. 2 where the 
proton chemical potential in pure neutron matter is shown for nuclear inter- 
action models from Fig.l. The neutron density profile n (r) is thus assumed 
to have a minimum at the center of the cell. 

The two other terms in eq.(4) describe the neutron background contri- 
butions to the energy. These represent the neutron Fermi sea energy and 
the curvature energy due to the gradient of the neutron distribution, re- 
spectively, in the Thomas-Fermi approximation. Here e(n(r)) is the local 
neutron matter energy per unit volume. The parameter is the curvature 
coefficient for pure neutron matter [2]. 

To decide which is the ground state configuration we compare the en- 
ergies Eq and El assuming the same number of neutrons in the cell. This 
means that the neutron density variation conserves the baryon number: 



In the next section the minimum of the energy difference AE = El — Eq is 
calculated in a simple variational approach and in Sect. 4 more sophisticated 
method is developed. 



We assume a simple trial form of the proton wave function and the 
neutron density variation. For the proton wave function we use a Gaussian 
form: 



Here Rp is the rms radius of the localized proton probability distribution. 
We treat this quantity as a variational parameter and minimize the energy 
difference AE with respect to Rp. 

Using the trial form of the proton wave function fyp (r) the energy dif- 
ference AE becomes 




(5) 



3. Simple estimate of the localized proton energy 





AE = 




(7) 
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The neutron density n (r) is chosen to be [2, 4]: 



n (r) = tin + a 



(8) 



Here a is the second variational parameter; a > corresponds to the neu- 
tron density enhancement around the proton and a < corresponds to the 
bubble in the neutron density near the proton. 

We calculate the energy difference AE, Eq.(7), for small proton fraction 
x, i. e. in the limit of large volume V. The first and the last terms were 
calculated assuming that the Wigner-Seitz cell radius Rc is much bigger 
than Rp, Rc 3> Rp. Denoting ^>* p (r) ^>p(r) = p (r) and expanding in y 
we have 



X 



p (r) ^fip (n N + ap (r) - a^j - hp (n N )j d 3 r 

/ p (r) {hp (n N + ap (r)) - fip (nN)) d 3 r- 
JV 



a vJy Pir 



dn 



{tin + ap (?")) d 3 r. 



(9) 



The integral in the last term does not depend on the cell volume so that 
this term vanishes in the limit V — > oo. Expanding in the same way the 
energy density, we obtain from the third term in Eq.(7) 



V 



e[n N + ap{r) - a—) - e {n N ) 



d 6 r = 



/ [e (n N + ap (r)) - e {n N )} d 3 r - a^N (njv) - 
Jv 

a T7 I {n N + ap (r)) - hn (njv)) d 3 r. 
V Jv 



(10) 



Here also the integral in the last term does not depend on the cell volume, 
since p (r) is a Gaussian, and this term vanishes for large V. The last term 
containing the coefficient of curvature Bn is easily evaluated to be: 
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The energy difference AE thus becomes: 
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{e (n (r)) - e (n N )} d 3 r - af i N (n N ) + ^ Qtt^J " ^B N a 2 . (12) 

We obtain physical parameters of the localized phase for a given neutron 
matter density tijy by a straightforward minimization of AE with respect 
to the two variational parameters a and Rp. 

The results of the calculations for the MS and FPR nuclear interactions 
are presented, respectively, in Figs. 3 and 4 where we show the energy dif- 
ference AE as a function of the proton distribution rms radius Rp for a 
few values of the neutron density. The curves are labeled with the value 
of the neutron matter density njy with subscript a. One can notice that 
for both MS and FPR interactions there appears a local minimum above 
a certain density, for the proton rms radius Rp in the range 1/m — 2fm. 
We have chosen the results for MS and FPR interactions only as examples 
of a general behaviour which is observed for all interactions we use in the 
calculations (more detailed account of our calculations will be presented 
elsewhere). With increasing neutron matter density un the depth of the 
minimum increases and above the threshold density the energy difference 
becomes negative. The negative value AE < means that the energy of 
the localized proton is lower than the energy of a nonlocalized proton and 
the localized proton state is preferred energetically. The behaviour of AE is 
very similar for all interactions we examine. This shows that the localization 
is not an effect of some specific interaction but rather is a general qualita- 
tive feature of the physical system we consider. Quantitative results, i. e. 
the localization density, the value of AE at the minimum and the localiza- 
tion radius Rp, depend on the specific interactions. The proton localizaton 
occurs at the lowest density for the Skyrme interactions, ni oc = 0.4/m -3 , 
and the energy difference AE displays the fastest decrease with the density. 
One can say that the localization is the strongest in this case. 

To understand better the localization mechanism it seems usefull to 
consider separately various contributions to the total energy difference. In 
Figs. 5 and 6 we show the proton contribution, Ep, to the energy difference 
AE, which consists of kinetic and potential terms. Here the minimum occurs 
at lower values of the proton rms radius Rp. One should keep in mind that 
the proton energy contribution represents a difference of the kinetic and 
potential energies of a localized proton and a single plane-wave proton. In 
the latter case the kinetic energy is zero. Next figures, Figs. 7 and 8, show the 
contribution of the neutron background to the total energy difference. This 
contribution is a monotonically decreasing function of the proton rms radius 
Rp. It grows very fast for low values of Rp. This rapid growth is similar to 
the behaviour of the gradient term contribution which is displayed in Figs. 9 
and 10. Thus a sum of these contributions also grows fast for low values 
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of Rp. Its values for a given radius Rp increase with the mean neutron 
density njy. 

The minimum of the total energy difference, AE, which is a sum of 
all contributions shown in Figs. 5-10, is thus a result of a delicate balance 
between repulsive contributions due to the neutron background, the proton 
kinetic term and the neutron curvature energy, and the attractive part of 
the proton interaction energy. Results of our calculations for a number of 
effective nuclear interactions show that such a minimum occurs in all cases 
above some density. One may thus conclude that the localization of proton 
impurity in the neutron matter is a general prediction of nuclear models of 
the class we consider here. 

The threshold density for proton localization, ni oc , depends also on the 
curvature coefficient B^ entering the gradient term and the proton effective 
mass mp, which are parameters of our model. In Fig. 11 we show how n/ oc 
changes with Bjy. One can notice that the localization density only weakly 
increases with the curvature coefficient in a wide range of its values. Also, 
the rms proton distribution radius at the threshold density R l p c increases 
slowly with increasing curvature coefficient Bn, Fig. 12. This fact is rather 
important in regard of validity of the Thomas-Fermi approach used in our 
model. 

The threshold density ni oc depends in a more sensitive way on the proton 
effective mass, mp, as shown in Fig. 13. For values of mp less than the 
bare proton mass the localization density increases. However, in the range 
QOOMeV — 938MeV which is most likely physically relevant to neutron stars, 
the threshold density changes by about 20%. 



Variational calculations of the localized proton energy presented in the 
previous section used the trial functions with only two variational param- 
eters, a and Rp. In this section we develop more advanced variational 
method which should give better estimate of the ground state energy of a 
localized proton. 

The energy difference AE, is a functional of two functions ^ (r) and 
n(r). The physical constraint is that the variation of the neutron back- 
ground conserves the baryon number. We should thus look for such func- 
tions ^(r) and n (r) that minimize the functional 



4. The self-consistent method 



/ [n (r) , vpp (r)] = AE-X 



Jv Uv 



f [n(r) -n N ]d 3 r-E f V* P (r) $> P (r) d 3 r - 1 , 



(13) 



8 



SZM PRINTED ON FEBRUARY 1, 2008 



where we explicitly include constraints of the baryon number conservation, 
Eq.(14), and the proton wave function normalization 



L 



y*p{r)^ P {r)d 6 r -1 = 0. (14) 

v 

The Euler-Lagrange equations corresponding to the functional (13) can 
be found easily. The differentiation with respect to gives the Schrodinger 
equation for the proton impurity: 



_ _L V 2 ^p (r) + [fi P (n (r)) - hp (n N )} V P (r) = E P ^ P (r) . (15) 
Imp 

Differentiation with respect to n (r) gives the second-order equation for the 
neutron density distribution n(r): 



S -^^n M *P (r) + W (n (r)) + 2B N ^1 - A _ 0. (16) 



The boundary conditions the functions \I/ (r) and n (r) obey at r — > oo are: 
n (r) = riAr, (r)| = 0. This allows us to identify the Lagrange multiplier 
A with the neutron chemical potential, 



A = (n N ) . (17) 

To calculate the cell energy we adopt for the proton wave function the 
Gaussian form used in the previous section and solve with this ansatz the 
equation (16). The rms radius of the proton probability distribution, Rp, 
is treated as a variational parameter. Numerical solutions of eq.(16) are 
presented in Figs. 14 and 15, where we show the neutron density distribu- 
tions n(r) obtained from equation (16). As one can notice the neutron 
background has somewhat different shape than that used in Sect. 3. In par- 
ticular, at higher mean neutron densities ni\r there appears a significant 
density enhancement at the well boundary which considerably strengthens 
the localization effect. With the simple method of Sect. 3 the neutron dis- 
tribution around the proton is a monotonically increasing function of the 
radius. 

Results of calculations of the energy difference AE with the self-consistent 
method are presented in Figs. 3 and 4 as curves labeled with the value of 
the neutron matter density only. In Figs. 5 - 10 the proton contribution and 
all components of the neutron background contribution to the full energy 
difference are shown together with those obtained in Sect. 3. The neutron 
background energy, Figs. 7 and 8, calculated with the self-consistent method 
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is well below the simple estimate of Sect. 3 for low values of the proton dis- 
tribution radius Rp. The reduction of the energy is even bigger for the 
gradient term contribution, Figs. 9 and 10. An opposite effect is observed 
for the proton energy contribution, Figs. 5 and 6, where the energy corre- 
sponding to the self-consistent method is higher than that calculated with 
simple trial functions in Sect. 3. 

5. Conclusions and implications 

Our self-consistent method gives lower energies of localized protons than 
the variational method used in sect. 3. In Figs. 3 and 4 we compare the energy 
difference AE obtained with both methods, for the MS and FPR models of 
nuclear interactions. As functions of the radius Rp the energy differences 
AE for the old variational method and for the selfconsistent method depart 
from one another only for small values of Rp. At Rp > lfm the curves 
for both methods in Figs. 3 and 4 are practically identical. However at 
small Rp the self-consistent method gives significantly lower energies. The 
minimum values of AE are considerably below those found with the old 
method and they occur at somewhat smaller radii. Our study indicates also 
that with the growing curvature coefficient the localization density grows. 
The results presented in Fig. 11 show that this growth is weak and the values 
corresponding to both methods are quite similar. Also, the dependence of 
the localization density on the proton effective mass, Fig. 13 is very similar 
for both methods. 

To conclude, the self-consistent calculations improve the estimate of the 
energy of the cell containing a localized proton, especially at small values 
of the rms proton radius Rp. The proton contribution Ep and the gradi- 
ent term contribution to AE are most affected by the new method. The 
ultimate goal is to calculate the proper wave function of the proton, which 
would give the true energy of the localized state. 

Results of our calculations for nuclear interactions we use indicate that 
the proton impurity in neutron star matter becomes localized at densities 
above 0.5 — 1.0/to~ 3 . The selfconsistent method gives lower energies of 
localized protons and smaller threshold localization densities than simple 
variational method with trial functions. This has important consequences 
for neutron stars as densities in this range correspond to inner core of neu- 
tron stars with masses exceeding one solar mass, M > 1M . In Fig. 16 we 
show neutron star masses corresponding to all nuclear interactions used in 
the calculations reported above. 
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Fig. 1. The proton fraction of the neutron star matter as a function of 
baryon number density for indicated nuclear interaction models. 
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Fig. 2. The proton chemical potential in pure neutron matter as a function 
of baryon number density for the same interactions as in Fig.l. 
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Fig. 3. The energy difference AE as a function of the proton rms radius 
for the Myers and Swiatecki interaction. The curves corresponding to the 
self-consistent calculations are labeled with the value of the neutron mat- 
ter density in [/m -3 ]. The curves labeled additionally with the letter a 
correspond to the simple method of Sect. 3. 
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Fig. 4. The same as in Fig.3 for the Friedman-Pandharipande-Ravenhall 
nuclear interaction. 
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Fig. 5. The proton contribution to AE for the Myers and Swiatecki model 
of nuclear interactions. Curves labeled as in Fig.3. 
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Fig. 7. The neutron background contribution to AE for the Myers and 
Swiatecki model of nuclear interactions. Curves labeled as in Fig.3 
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Fig. 8. The same as in Fig. 7 for the Friedman-Pandharipande-Ravenhall 
interaction. 
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Fig. 9. The gradient term contribution to AE for the Myers and Swiatecki 
nuclear interactions. Curves labeled as in Fig.3 
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Fig. 10. The same as in Fig.9 for the Friedman-Pandharipande-Ravenhall 
interactions. 
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Fig. 11. The threshold density for proton localization versus the curvature 
coefficient B^. The curves MSI and FPR1 correspond to the simple method 
of Sect. 3. The curves MS2 and FPR2 correspond to the self-consistent 
calculations. 
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Fig. 12. The rms proton distribution radius at the threshold density as a 
function of the curvature coefficient B^. Curves labeled as in Fig.ll 




Fig. 13. The threshold density for proton localization as a function of the 
proton effective mass. Curves labeled as in Fig.ll 
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Fig. 14. The neutron density distribution obtained from Eq.(16) for in- 
dicated neutron matter densities (in [fm^ 3 ]) for the Myers and Swiatecki 
interactions. The localized proton distribution np is also shown. 




Fig. 15. The same as in Fig. 14 for the Friedman-Pandharipande-Ravenhall 
interactions. 
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Fig. 16. Neutron star masses (in solar units) as functions of the central 
density for nuclear interaction models from Fig.l. 




